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天 中 平原 PM2s 浓 度 周期 性 变化 规律 及 成 因 分 析 
Mt 


摘 要 : 采用 Morlet 小 波 分 析 法 和 交叉 小 波 分 析 法 ,对 关 


， ERR, wh 哺 
(中 国人 民 大 学 环境 学 院 , 北 京 


100872 ) 


Pp 五 市 2013 年 11 月 15 日 至 2019 年 3 月 15 日 期 间 的 PM; 


浓度 周期 性 变化 特征 及 成 因 进行 分 析 。 结 果 表 明 :(1) 关中 五 市 的 PM; 浓 度 变 化 具有 显著 区 域 性 和 周期 性 , PM; 


浓度 在 夏季 呈现 低 浓 度 稳 定 状 态 , 秋 冬季 呈现 高 浓度 状态 ， 


呈现 出 365 d 左 右 、.36 d 左 右 以 及 14 dd 左右 的 周期 性 变 


化 规律 .(2) 秋冬 季 , 关 中 五 市 的 PM;; 浓 度 以 14 qd 左右 的 周期 性 变化 为 主 。PM;; 浓 度 14 d 左 右 的 周期 性 变化 与 风 


向 、 风 速 和 气压 三 项 气象 因子 基本 为 负 相 关 关 系 , 与 湿度 呈现 落后 1/8 个 相位 周 共 


,与 气温 基本 无 相关 关系 ,(3) 关 


中 五 市 秋冬 季 PM; 浓 度 周期 性 变化 主要 是 由 运城 三门峡、 洛阳 和 关中 城市 等 区 域 形成 的 局 地 静 稳 气流 (该 类 气流 
占 21.08%, 在 该 类 气流 下 PM; 不 断 积 累 ) 和 从 蒙古 国 、 内 蒙古 等 地 输入 的 气流 (该 类 气流 占 23.86% ,在 该 类 气流 下 


积累 的 PM2s 得 到 清除 ) 周 期 性 交替 所 致 。 


关键 词 : 关中 平原 ; PM;; 浓 度 ; 周期 性 变化 ; 成 因 分 析 


关中 平原 位 于 陕西 省 中 部 地 区 ,位 于 渭南 、 西 
安 .宝鸡 .咸阳 和 铜川 五 市 (也 称 关中 五 市 ) ,被 北山 
山系 、 和 秦岭 以 及 晴 山 等 高 海拔 山脉 环绕 ,形成 了 独 
特 的 地 形 气候 ” 。 随 着 关中 五 市 经 济 的 迅速 发 展 ， 
人 为 活动 源 排放 污染 物 不 断 增 加 ,秋冬 季 重 污染 天 
气 频 发 “ 。 针 对 关中 平原 的 PM:: 污 染 问题 ,多 数 研 
究 集中 在 应 用 WRF 等 污染 扩散 模型 对 成 因 进行 分 
析 m" ,从 季 闻 、 月 等 长 时 间 尺 度 变换 的 角度 对 PM 
变化 规律 的 研究 较 少 。 分 析 PM;; 周 期 性 变化 趋 
势 , 有 利于 优化 PMs 预 测 方法 和 更 好 地 治理 PM 
污染 ”。 小 波 分 析 是 研究 时 间 序 列 中 周期 性 现象 的 
有 效 方法 ,尤其 在 识别 随时 间 变 化 的 潜 频 率 方面 效 
FYE. H 1980 年 代 初 以 来 ,小 波 分 析 大 量 应 用 
于 信号 和 图 像 处 理气 象 水 文 变 化 地球 物 理 和 经 
济 学 等 方面 ""。 最 近 几 年 , 李 梓 铭 '“、 王 海 鹏 "、 
周 静 ”等 采用 小 波 分 析 方 法 分 别 对 北京 .兰州 、 太 
原 地 区 的 PM; 长 时 间 的 浓度 变化 趋势 进行 分 析 , 识 
别 PM;;s 变 化 规律 。 本 人 研究 采用 趋势 分 析 、 小 波 分 
析 、 交 又 小 波 分 析 等 方法 ,人 研究 关中 平原 的 PM;; 发 
展 规律 .变化 特征 、 影 响 因 素 ,将 有 助 于 了 解 关 中 平 
原 PM;; 污 染 分 布 特征 、 变 化 规律 及 成 因 。 


收 稿 日 期 : 2020-04-03; 修订 日 期 : 2020-06-04 


1 研究 方法 


1.1 数据 来 源 

PM:: 数 据 和 风速 ` 风向 气温 .气压 ,湿度 等 气 
象 数据 来 自 关 中 五 市 空气 质量 国 控 监 测 点 。 研 究 
时 段 为 2013 年 11 月 15 日 至 2019 年 3 月 15 日 ,气象 
数据 从 2014 年 11 月 15 日 至 2019 年 3 月 15 日 , 复 盖 
关中 五 市 6 个 冬 防 期 (11 月 15 日 至 次 年 3 月 15 
日 )。 后 向 轨迹 模式 (HYSPLIT) 聚 类 分 析 所 用 的 气 
象 场 数 据 来 自 美国 国家 海洋 和 大 气管 理 局 (NOAA) 
(ftp://arlftp.arlhq.noaa.gov/pub/archives/reanalysis/ ) 。 
使 用 R 软件 WaveletComp 包 进行 Morlet 小 波 分 析 和 
交叉 小 波 分 析 ”。 
1.2 Morlet 小 波 分 析 方 法 

经 典 的 傅立叶 变换 可 以 把 时 间 序 列 变换 到 频 
域 ,但 只 能 单独 从 时 域 或 频 域 分 析 时 间 序 列 。 小 波 
变换 克服 了 时 域 和 频 域 不 能 兼顾 的 缺点 ,可 以 从 时 
域 和 频 域 两 个 角度 同时 分 析 时 间 序 列 。Morlet 小 波 
通过 频 移 保 持 形 状 , 在 时 间 分 辩 率 基本 无 损失 的 情 
况 下 ,对 不 同 频带 的 波 进行 有 效 分 离 和 级 数 重 构 ， 
有 具有 良好 的 时 间 聚 集 性 、 较 高 的 频率 分 辨 率 、 包 含 
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相位 信息 等 特点 ”中 。 本 研究 选用 Morlet 小 波 对 关 
中 平原 的 PM 的 周期 性 进行 分 析 。Morlet 小 波 的 时 
域 , 频 域 表 达 式 为 : 
h(t) = ae ee ( 1 ) 

式 中 :1 为 时 间 ; o, 是 无 量 纲 频率 。 当 wo =6, 小 波 尺 
度 * 与 传 里 时 周期 基本 相等 ,尺度 项 与 周期 项 可 以 
HERR 

将 小 波 功率 谱 |W,(9| 在 某 一 周期 上 进行 时 间 
平均 ,可 以 得 到 一 幅 展 示 时 间 系 列 在 某 一 尺度 上 的 
波动 特征 及 其 随时 间 变 化 的 小 波 全 谱 , 小 波 全 谱 能 
表明 时 间 系 列 真实 功率 谱 的 无 偏 一致 估计 "oo2n。 
由 于 小 波 全 谱 可 以 显示 出 背景 谱 ,通过 小 波 全 谱 
图 ,可 以 清晰 辨别 时 间 系列 的 周期 波动 特征 和 
强度 。 


W (s)= 15) LAOI (2) 


RP: n 为 总 时 间 个 数 ; WO) 为 连续 变换 小 波 函数 。 
1.3 交叉 小 波 分 析 

交叉 小 波 分 析 可 以 比较 两 个 时 间 序 列 多 和 的 
频率 ,得 出 两 个 序列 在 某 些 时 间 段 内 的 共振 周期 和 
相位 情况 。 两 个 时 间 序 列 的 交叉 小 波谱 图 XWT 为 
WY=W*W”, 交 叉 小 波 功率 为 |W”| 。 同 样 ,交叉 
小 波 功率 所 反映 的 周期 随时 间 变 化 的 性 质 ,可 以 用 
全 谱 图 进行 可 视 化 展示 "mi 。 交 又 小 波 宕 的 置信 度 
可 以 从 两 个 卡 方 分 布 的 乘积 的 平方 根 得 出 ,其 置信 
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, DLAKO) 


sw fx sw | 

交叉 小 波 变 换 包 含有 关 两 个 序列 的 同步 性 信 

息 ,这 些 信息 取决 于 一 个 序列 的 任何 周期 分 量 相 对 

于 男 一 个 序列 相应 分 量 的 局 部 相位 差 情 况 。 相 位 

差 可 以 转换 为 间隔 [-T, T] 中 的 某 个 角度 的 差 值 , 相 

位 差 的 大 小 可 以 反映 两 序列 的 济 后 特性 , 正 负 可 以 

反映 两 序列 的 主导 系列 ”。 相 位 差 公 式 如 下 : 
人 


R[W,""(s)| 


Px = ws (4) 


2 结果 与 讨论 


21 关中 五 市 PMs 浓 度 变化 趋势 

根据 关中 五 市 PMs; 浓 度 从 2013 年 11 月 15 日 至 
2019 年 3 月 15 日 的 月 尺度 浓度 变化 趋势 (图 1)， 
PM; 浓 度 存 在 明显 的 季节 性 波动 , 冬 防 期 为 当年 浓 
度 变 化 的 波峰 期 ,夏季 为 浓度 波 谷 期 。 在 Ps 浓度 
的 波动 过 程 中 ,风速 、 风 向 等 气象 因素 对 于 PM; 的 
生成 .转移 和 去 除 有 密切 的 关系 站 。 根 据 图 2, 关 中 
五 市 PW 大 于 80 mg m3 时 ,风向 分 布 基 本 在 东北 
风 至 北 风 之 间 ,风速 分 布 在 0~4m:'s 2], 1 PM25 
浓度 值 超过 110 mg-:m ;时 ,基本 处 于 微风 静 稳 天 
气 。 由 于 关中 五 市 PY; 浓度 波动 存在 一 致 的 变化 
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图 1 关中 五 市 PW; 浓度 变化 趋势 


Fig. 1 Trend of PM>; concentration in five cities of Guanzhong 
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图 2 关中 五 H PM;; 浓 度 与 风速 风向 分 布 


Fig2 PM;s concentration and wind speed and direction in five cities of Guanzhong 


趋势 ,高 浓度 PM;; 发 生 时 风速 、 风 向 等 气象 条 件 基 
本 一 致 ,推断 关中 平原 的 PM;; 浓 度 变化 存在 较 强 区 
域 性 。 

2.2 关中 五 市 PM2: 周 期 变化 规律 

为 分 析 2013—2019 年 关中 五 市 PM:s 的 时 间 周 
期 特征 ,分 别 对 关中 五 市 2013 一 2019 年 日 均 PM:: 浓 
度 的 时 间 序 列 作 Morlet 小 波 变换 ,进行 小 波 功率 谱 
图 分 析 ( 图 3)。 在 小 波谱 图 中 ,白色 阴影 包围 的 蓝 
色 区 域 为 小 波 边界 效应 影响 锥 ,其 周期 性 特征 存在 
较 高 不 确定 性 ,白色 粗 实 线 所 包围 的 红色 区 域 表示 
通过 显著 性 水 平 为 a=0.1 的 噪声 检验 ,其 对 应 的 时 
间 段 具有 显著 周期 性 。 

从 小 波谱 图 中 发 现 , 关 中 五 市 Ps 浓度 在 冬 防 
期 间 均 存在 4~32 d 的 显著 周期 性 变化 ,与 关中 平 
原 在 冬 防 期 间 周 期 性 出 现 雾 儿 天气 的 现象 吻合 。 
结合 小 波 功率 曲线 图 可 以 看 出 ,关中 五 市 PY; 存在 
两 个 显著 的 周期 ,第 一 主 周 期 均 为 365 d 左 右 , 第 二 
主 周期 均 为 14 d 左 右 ,这 两 个 周期 均 通过 了 90% 置 
信和 度 检 验 。 同 时 ,咸阳 市 西安 市 和 渭南 市 36 d 左 
右 , 宝 鸡 市 和 铜川 市 32 4 左右 的 周期 变化 也 通过 了 
90% 的 置信 度 检验 。 此 外 ,咸阳 市 .西安 市 和 渭南 市 
存在 6d 左 右 的 周期 变化 ,宝鸡 市 和 铜川 市 存在 7d 


节 变 化 以 及 短期 变化 特征 。 以 西安 市 为 例 , 绘 制 第 
一 主 周期 365 d、 第 二 主 周期 14 d、 季 市 性 周期 36 d 
和 星期 效应 7d 时 间 尺 度 下 的 小 波 时 序 图 ,PM;s 浓 
度 在 夏季 保持 低 浓 度 稳定 状态 ,春季 有 小 范围 波 
动 , 秋 冬季 呈现 出 强烈 的 震动 性 ,从 36 d 到 7d 之 
间 , 震 荡 周 期 越 来 越 短 , 频 率 越 来 越 高 ,在 12 月 至 次 
年 1 月 期 间 达 到 最 大 振幅 。 从 图 4 可 知 , 在 2016 一 
2017 年 冬 防 期 间 ,7d 和 14 d 周 期 下 的 PM;; 浓 度 波 
动 幅度 明显 高 于 其 他 年 份 。 
2.3 PM 浓度 周期 性 变化 与 气象 因素 之 间 的 关系 
通过 以 上 分 析 可 知 ,关中 五 市 的 PW; 浓度 周期 
变化 呈现 区 域 性 特征 。 以 西安 市 为 例 ,分 析 风 向 、 
风速 .气压 、 气 温和 湿度 对 PM; 周 期 性 的 影响 。 将 
西安 市 PW: 浓度 的 时 间 序 列 分 别 与 各 气象 因子 作 
交叉 Morlet 小 波 变换 ,进行 交叉 小 波 功 率 谱 图 分 析 ， 
进一步 分 析 PM;; 浓 度 周期 。 在 交叉 小 波谱 图 中 (图 
5) ,白色 阴影 包含 的 蓝 色 区 域 为 小 波 边界 效应 影响 
锥 ,白色 粗 实 线 包围 的 红色 区 域 表示 通 过 显著 性 水 
平 为 a=0.1 检 验 , 箭 头 表 示 两 序列 之 间 的 相位 关系 ， 
“一 "表示 PM; 与 气象 数据 之 间 为 同 相 位 ,二 者 为 正 
相关 关系 ,“<*-” 表 示 PM; 与 气象 数据 之 间 为 反 相 
位 ,二 者 为 负 相 关 关 系 ,“ ”和 “1 ”分 别 表示 PMs 


左右 的 周期 变化 ,关中 平原 的 PM,; 变 化 整体 存在 星 
期 效应 。 
分 析 关 中 平原 2013 一 2019 年 PM., ,长 期 变化 E 


与 气象 数据 相 比 提前 和 落后 1/4 个 周期 , 非 线 性 
相关 。 
从 图 5 可 以 看 出 , 冬 防 期 间 PM2:: 浓 度 与 风向 、 风 
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图 3 关中 五 市 PMzs 小 波谱 图 及 小 波 功 率 曲线 分 布 


Fig. 3 PM,s Wavelet spectrum and wavelet power curve of five cities in Guanzhong 
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注 : 图 中 灰 线 为 实际 PM2: 浓 度 , 黑 线 为 特定 时 间 尺 度 下 小 波 系 数 
4 特定 时 间 尺 度 的 小 波 时 序 


Fig. 4 Wavelet time series of specific time scale 
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图 5 西安 市 PM2s 浓 度 与 气象 因子 的 交叉 小 波谱 


Fig. 5 Cross wavelet spectrum of PM2s concentration and meteorological factors in Xi’ an 


速 .气压 湿度 在 周期 8~32 d P-L a tee , 形 
成 了 14 d 左 右 的 第 二 主 震 荡 周 期 .7 d 左 右 的 第 三 震 
荡 周 期 36 4 左右 的 震荡 周期 。PM:: 与 风向 、 风速、 
气压 基本 为 负 相关 , 与 湿度 呈现 落后 1/8 个 相位 周 
期 。PM :浓度 与 风向 ` 气压 湿度 在 2016 一 2017 年 
冬 防 期 间 , 在 7d 左 右 和 14 d 左 右 的 共振 程度 明显 
高 于 其 他 年 份 , 说 明 关 中 平原 2016 一 2017 年 冬 防 期 
间 的 PMzs 浓 度 较 其 他 年 份 高 ,受到 气象 因素 的 不 利 
影响 。 
2.4 周期 性 变化 成 因 分 析 

在 北半球 冬 半年 平流 层 中 ,存在 显著 的 10~ 20 
d 双 震荡 周期 中 和 30~60d 的 季节 震荡 周期 ”等 
低频 振荡 ,这 些 震荡 导致 了 风速 .风向 .气压 等 气象 
要 素 的 变化 ,在 水 平方 向 上 由 东 向 西 传播 "。 通 过 
后 向 轨迹 聚 类 分 析 , 可 以 将 具有 相似 气象 因素 的 气 
流 进行 聚 类 分 组 ,研究 相似 地 理 起 源 的 数据 所 对 应 


的 污染 物 传输 情况 ”2 。 通 过 以 西安 市 为 起 点 (图 
6), 对 2013 年 11 月 15 日 至 2019 年 3 月 15 日 期 间 的 
近 地 面 (10 m)5 d 时 长 的 后 向 轨迹 气流 以 季节 尺度 
(春季 :3 一 5 月 ,夏季 :6 一 8 月 ,秋季 :9 一 11 月 ,冬季 
12 一 次 年 2 月 ) 进 行 聚 类 ,发现 西安 市 四 季 均 存在 从 
运城 三门峡 洛阳 、 渭 南 等 地 输送 的 西风 西北 风 
局 地 气流 ,以 及 从 天 水 .陇南 和 宝鸡 .成 阳 等 地 输送 
的 东风 、 东 北 风 长 距离 气流 。 在 冬季 ,来 自 运 城 .三 
门 峡 、 洛 阳 等 由 东 向 西 传播 的 C4 类 气流 占 关中 平原 
冬季 气流 的 21.08%( 以 西安 咸阳、 渭南 、 宝 鸡 和 铜 
川 为 起 点 ,该 类 气流 分 别 占 冬 季 气 流 的 22.7%、 
23.3% 、22.3% 、17.7% 和 19.4%)。 对 关中 平原 冬季 
聚 类 后 的 C1.C2、.C3、C4 四 类 后 向 轨迹 气流 对 应 
PM; 浓 度 进 行 分 析 ( 图 7), 当 C4 类 局 地 气流 出 现 
时 ,整个 关中 平原 PM:: 呈 现 高 浓度 水 平 ,结合 图 2 可 
知 , 高 浓度 PM; 出 现时 关中 平原 气象 条 件 处 于 静 稳 
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图 6 西安 市 不 同 季节 下 的 近 地 面 大 气 输送 通道 
Fig. 6 Near surface atmospheric transport channels in different seasons in Xian 
[ 一 Cl 类 气流 。 ”分 别 占 冬 季 气 流 的 25.5% 21.8% , 25.5% 、24.7% 和 
200 上 一 0 一 C2 类 气流 E E ae 
| 到 C3 类 气流 21.8% ) ,与 醉 文博 和 黄 光 球 等 ”的 研究 结论 基本 
= < CAR 吻合 。 同 时 从 图 7 可 以 看 出 ,在 冬季 各 类 气象 条 件 
% 160} F, PEZAR PIT PM AS HH AT DX SIS YE, 
= 140} PM;; 浓 度 均 高 于 关中 其 他 地 市 ,也 在 一 定 程度 上 说 
总 1204 明了 西安 和 咸阳 两 市 两 地 城市 PMss 本 地 源 排放 强 
2 rool 度 相对 排放 量 要 高 于 关中 其 他 地 市 。 
综 上 分 析 , 可 以 看 出 由 于 关中 平原 地 形 特殊 ， 
60l Be 结合 气象 影响 中 和 污染 物 传输 及 排放 ,使 得 关中 平 
en on ae a 原 的 PM; 浓 度 在 冬季 形成 了 14 d 左 右 的 双 震 荡 周 
关中 五 市 期 以 及 36 d 左 右 的 季节 震荡 周期 。 同 时 , 闭 加 本 地 污 
图 7 关中 五 市 冬季 不 同 后 向 轨迹 聚 类 下 的 PY; 


浓度 分 布 
Fig.7 PM;s concentration under different backward trajectory 


clusters in five cities of Guanzhong in winter 


状态 ,可 以 推测 在 此 类 不 利 气象 条 件 下 ,本 地 排放 
量 较 大 ,超过 相应 气象 条 件 下 的 环境 容量 ,使 得 
PM:: 浓 度 升 高 。 

当 冬 季 出 现 途 经 蒙古 国 内蒙、 宁夏 等 地 的 长 
距离 输送 C2 类 气流 时 ,该 气流 对 关中 平原 的 PM; 
进行 清除 ,使 得 整个 关中 平原 的 PM; 浓 度 在 较 低 水 
平 (图 7) ,该 类 气流 占 关中 平原 冬季 气流 的 23.86% 
(以 西安 咸阳、 渭南 .宝鸡 和 铜川 为 起 点 ,该 类 气流 


染 源 的 排放 影响 后 ,形成 了 7 d 左 右 的 星期 效应 Pi。 
3 结论 


(1) 关中 平原 五 个 城市 PM;; 浓 度 变化 规律 基本 
保持 一 致 ,具有 显著 区 域 性 和 周期 性 。 夏 季 保 持 稳 
定 低 浓度 状态 ,秋冬 季 呈 现 季 节 性 高 浓度 状 周期 性 
变化 。PM; 浓 度 在 夏季 呈现 低 浓 度 稳定 状态 ,秋冬 
季 呈 现 高 浓度 状态 。 其 中 ,365 d 左 右 、36 d 左 右 以 
及 14d 左 右 的 主 震 荡 周 期 通过 了 90% 的 置信 度 
检验 。 


(2) 秋冬 季 , 关 中 五 市 的 PM; 浓 度 以 14 dA 
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的 周期 性 变化 为 主 。14 d 左 右 的 PM:: 浓 度 周 期 性 China Environmental Science, 2019, 39(8): 3539-3548. | 
7S fe 3 BE SE RU RL] FR A BE DS TTB PN, BIR, SP PDE PM AAEM DORR 
进行 驱动 。 结 合 后 向 轨迹 聚 类 分 析 , 得 出 其 原因 主 [J]. 生态 学 杂志 , 2019, 38(8): ane ea pe 2i Li Yungi, 
z RRE he 5 ‘ 3 ais Zhao Yonghua, et al. Ground concentration distribution and source 
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Analysis on the periodic variation of PM., concentration and its 
causes in Guanzhong Plain 


JIA Ce, WANG Ya-qiong, YANG Xiao 
(School of Environment & Nature Resources, Renmin University of China, Beijing 100872, China) 


Abstract: To analyze the characteristics and causes of PM2s concentration changes in five Guanzhong cities from 
November 15, 2013 to March 15, 2019, Morlet wavelet analysis and cross wavelet analysis were used. The results 
show that (1) the changes in PM2s concentration in five Guanzhong cities were markedly regional and periodic. The 
PM2s concentration was low in summer and high in autumn and winter, with periodic changes every 365 d, 36 d, and 
14 d, respectively. (2) In autumn and winter, the PM2s concentration in five Guanzhong cities changes periodically 
every 14 d. The cyclical change in PM;s concentration every 14 d was negatively correlated with the three meteorological 
factors of wind direction, wind speed, and air pressure, and was lagging behind by 1/8 phase period with humidity, 
with no correlation with temperature. (3) The periodic change in PM2s concentration in five Guanzhong cities in 
autumn and winter was mainly due to the local static airflow formed by Yuncheng, Sanmenxia, Luoyang, and 
Guanzhong cities (this type of airflow accounts for 21.08% , and under this type of airflow, PM25 continuously 
accumulates) and the airflow imported from Mongolia, Inner Mongolia, and other places (this type of airflow 
accounts for 23.86%, and the PM2s accumulated under this type of airflow is removed) periodically. 


Keywords: Guanzhong Plain; PM2s concentration; periodic change; cause analysis 


